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Abstract 

Two- parameter sets of solutions to the classical field equations in the mass- 
less (p 4 model and SU(2) gauge theory are found, each solution presumably 
describing a multi-particle instanton-like transition at high energy. In the limit 
of small number of initial particles, the probability of the transition is sup- 
pressed by exp(— 2 So), where So is the instanton action. 
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1 Introduction 



For quite a while, there has been a considerable interest in the problem of non- 
perturbative scattering at parametrically high energies in weakly coupled theories. 
In the context of the standard model, this problem is usually associated with the 
anomalous baryon number violation which might be observable in particle collisions 
in 10-100 TeV region. Despite the considerable effort made since the first quantitative 
attempts to estimate the corresponding cross section |l], Q , the complete understand- 
ing of this problem has not yet been achieved (for a review of the present status see 
refs.[|,l). Although the total cross section at high energies has been strongly argued 
[HI to have the semiclassical form 



ffW = exp{^(|)} 



r 

where g is the coupling constant and E ~ My/ / g 2 is the characteristic energy scale 
(the sphaleron scale), the explicit procedure for evaluating the function F(E/Eq) is 
still unknown, the main difficulty consisting in the semiclassical treatment of high 
energy initial particles. 

To overcome this difficulty, there has been recently proposed an implicit approach 
|], [7] which makes use of the idea to approximate the two-particle cross section in 
the one-instanton sector by the probability of a multiparticle instanton-like transition 
for a small number of initial particles. If the number of initial particles scales like 
n = v/g 2 , where v is a (small) numerical constant, the latter probability is explicitly 
semiclassical, 

a(E,n) = expl±-F(^,is)\ , (2) 



while the function F(E/E , u) is calculable by semiclassical methods || [7|, § for any 
particular coherent initial state containing n = vjg 2 particles. 

There are two possible ways of extracting information about the two-particle cross 
section from this sort of considerations. First, one may study the probability of the 
transition summed over (or, equivalently, maximized with respect to) some subset 
of fixed-energy initial states. The obvious, but by no means unique choice could be 
the set of all n-particle states of fixed energy. This maximum probability is again 
semiclassical, i.e. has the form (eq.(|[)), where the exponent depends, of course, on 
the particular choice of the subset of the initial states. By construction, for all sets 
that include the particular state where two particles carry almost all the energy, the 
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corresponding function F(E/Eq, v) sets an upper bound for the logarithm of the two- 
particle cross section, F(E/E ). The advantage of this way is that it allows not to 
deal with (probably irrelevant) details of the initial state. 

An alternative way is to consider the transition from some particular coherent state 
| a) containing n = vjg 2 particles. The perturbative calculations in three first orders 
in the electroweak theory indicate that, in the limit v — > 0, the corresponding function 
F a (E/E , v) is smooth and independent of the details of the initial state ||. If this 
property held in higher orders, the limiting value F a {E j Eq^v — > 0) would coincide 
with the small v limit of the maximum probability discussed above and, furthermore, 
would coincide with the logarithm of the two-particle cross section, F(E/E Q ). 

Clearly, the first way is preferable as it gives the upper bound on the two-particle 
cross section regardless of whether the assumption about the independence of the 
probability of the precise form of the initial state is justified. Unfortunately, it requires 
solving the field equations with particular boundary conditions f8|, which is almost 
certainly beyond the scope of analytical methods. The more realistic way could be to 
find any solutions to the field equations and then check if they correspond to some 
instanton-like transition with small number of initial particles. This is our strategy 
in the present paper. 

There are two requirements which necessarily have to be satisfied by the solution 
saturating the total probability of instanton-like transition from some initial state ||: 

i) The solution must be instanton-like, i.e. must interpolate between two different 
asymptotic regions separated by a potential barrier. The underbarrier part of 
the field configuration is described by the evolution in the imaginary time (the 
"duration" of the Euclidean stage T is infinite for the instanton and finite for 
finite energy transitions). So, the solution we are interested in should exist on 
the contour ABCD of fig.l in the complex time plane and satisfy the complex- 
ified field equations along this contour. Since the initial and final asymptotic 
regions, A and D in fig.l, must lie on different sides of the barrier, there must 
be singularities between the line AB and the negative part of the real time 
axis, which prevent from doing trivial analytic continuation between these two 
regions. 

ii) The solution must be real on the real time axis. This requirement is valid as 
long as the probability summed over the final states is considered. Note that 
it can be viewed as one of the two boundary conditions necessary to specify a 
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unique solution to the field equations. In what follows we will impose a stronger 
condition, namely requirement that the solution has a turning point at t = 0, 
i.e. is also real at the imaginary time axis. 

Any solution satisfying these two conditions corresponds, in general, to the instanton- 
like transition from some initial state to the state on the other side of the barrier 
which maximizes the transition probability. In other words, the solutions of this type 
saturate the total probability 

a a (E)~J2\(b\SP E \a)\ 2 (3) 

b 

at some initial state \a). Here Pe is the projector onto fixed energy E. 

Given the solution satisfying the conditions i) and ii), one can easily find the 
characteristics of the instanton-like transition it corresponds to. Assuming that fields 
are free in the asymptotic regions A and D, one writes 

&(M) = i^(/ k e-**i + /_ k e«*"), (4) 

cf> f (k,t) = ^^(6 k e-^ + 6*_ k e^), (5) 

where r\ = t — iT is a real parameter along the line AB. Note that because of the 
condition ii) the amplitudes &k and 6 k are complex conjugate to each other, while 
/k and /k are not. The most probable final state of the transition can be read off 
from eq.(f|). This is the coherent state defined by the complex variables 6k, \b) = 
exp(/<ifc&k6k)|0). Correspondingly, the number of particles in the final state and its 
energy are 

n f = Jdkblb k , 

r (6) 

E = dku k blbk . 

The complex variables characterizing the initial coherent state \a) can be obtained 
from the amplitudes /k entering eq.(f|) as follows, 

/k = a k e^ AT . (7) 

The presence of the additional parameter AT reflects the fact that the coherent 
states | Ok) and P^l^k) differ only by normalization |H|. The number of particles in 
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the initial state depends only on the difference T — AT, 

n % = J dfeajok = J dfc|/ k |V 2 ^ AT . (8) 

The latter difference can be found from the energy conservation, 

J dkco k b k b* k = J dku^ai = J dkco k \f k \ 2 e- 2uJkAT . (9) 

Finally, the probability of the transition is given by the formula 

a = exp (2E(T - AT) - 21m S((/>) + J dfc(/k/k - a*a k )| , (10) 

where S((j)) is the action calculated along the contour ABCD of fig.l. 

The examples of the solutions satisfying the requirements i) and ii) were recently 
found in the two-dimensional 0(3) a- model In this paper we generalize the 

method of ref . |P_0|| to the four- dimensional case, namely, to the massless 4 model 
and the massless SU(2) gauge theory. The latter can be viewed as the high energy 
limit of the electroweak theory. In sections 2 and 3 devoted to the 4 model and 
SU(2) theory, respectively, we show that in both cases there exists a two-parameter 
family of the solutions satisfying the requirements i) and ii) . One of these parameters 
is a trivial scale parameter which is due to the conformal symmetry of the models. 
Another parameter, e, is dimensionless. At small values of this parameter the corre- 
sponding instanton-like transition is characterized by the number of initial particles 
much smaller than the number of final particles, both being much less than 1/g 2 , the 
characteristic number of particles in the sphaleron solution. The probability of the 
transition calculated by means of eq. (ffD) is exp(— 2S inst + 0(e 2 )), i.e. is suppressed 
by the same factor as at zero energies. The existence of such solutions may indicate 
that the probability of the instanton-like transitions from the initial states contain- 
ing small number of particles is suppressed at high energies by the same suppression 
factor as at zero energy. Further discussion of this possibility is contained in sect. 4, 
where we also present our conclusions. 



2 </> 4 model 

As a simplest four-dimensional model, consider the 4 theory with negative coupling 
constant. This model was suggested as a suitable testing ground for studying in- 



stanton transition at high energies [11], [12]. The Euclidean Lagrangian of the model 
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reads 



i 



(11) 



where A is positive. In this model, the state = is metastable with respect to 
non-homogeneous fluctuations. The decay of this state is described by the instanton 
solution 



of the instanton solution 4>i ns t- These processes share many common features with 
instanton-like transitions in the electroweak theory (such as exponential growth of 
the total cross section at small energies) and formally may be used as a laboratory 
for studying the instanton effects in particle collisions. 

As discussed in the Introduction, the probability of the instanton-like transition is 
saturated by the classical solutions satisfying the requirements i) and ii). To find such 
solutions one can use the method of ref.||10| based on the conformal invariance. First, 
one looks for an 0(4) symmetric solution to the Euclidean field equations. Making 
use of the spherically symmetric ansatz 



where £ = ln(y/R) and R is a parameter with the dimension of length, we find the 
equation for /(£) which is formally the classical equation of motion for a particle 
moving in one-dimensional double-well potential, 




(12) 



The action of the instanton is 





(13) 



d 2 f 



f ~ 2/ 3 . 



The solution to this equation with the turning point at £ = reads 




(14) 



where k = (1 — e 2 ) 1 / 2 and e is a remaining free parameter which varies in the range 
< e < 1. In eq.(|n|) K(k) is the complete elliptic integral. The instanton solution 
fll2"|) is reproduced at e = (k = 1). 
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In order to obtain the solution to the field equations which has the turning point 
at t = 0, we perform the conformal transformation to the new variables x^, 



= (y + a y (Vi* + q m) ~ a M » ( 15 ) 



with a M = (i?, 0,0,0), which maps the sphere y 2 = R 2 onto the upper half of the 
x-space. The solution (j!3|) becomes 

(X) - 



J(( a;o _ j R)2 + r 2)(( a . + i2 )2 + r 2) 



where r = |«c|. It is convenient to perform the analytical continuation to the Min- 
kowskian domain in this equation. Finally we obtain 

m = Mm (16) 

V\ yj((t - iR) 2 - r 2 )((t + iRY - r 2 ) ' 

where 

1^ {t-iRf-r 2 



2 (t + iR) 2 - r 



2 



Due to the conformal invariance of the 4 model, the field configuration (16) is again a 
solution to the field equations. Furthermore, since /(£) has the turning point at £ = 0, 
the field flOp is clearly symmetric under the time reversal t — > — t and, therefore, has 
a turning point at t — 0. 

Consider now the structure of the solution ([16]) in the complex time plane. As 



can be seen from eq. ( 16|) , for each r there exists an infinite set of singularities located 
as shown in fig. 2. At r — > oo and small e the positions of these singularities can be 
found analytically, 



2(2n+l) N 



t n = -r + iR\l-2{-) j . (17) 

It can also be shown that there exists a contour Imt = T such that at all r it passes 
above the first singularity in this series and below the next one. Since our solution 
is real on the real time axis and has a turning point, both requirements i) and ii) are 
satisfied. 

Now we have to determine what kind of scattering process the above solution 
corresponds to. For this purpose we consider the asymptotics of the field fll6|) in the 
regions t — > %T — oo and t — > oo, which correspond to the initial and final states, 
respectively. In the asymptotic regions the solution becomes a linear superposition of 
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the plane waves, the coefficients in this superposition determining the initial coherent 
state and the most probable final state of the process. Making use of the explicit 
solution ([TED, it is easy to check that in the asymptotic regions the field is concentrated 
at r — \t\ ~ 1: (j) falls rapidly as one moves off the light cone and becomes negligibly 
small when r — \t\ ^> 1. The magnitude of the field on the light cone decreases like 
1/t for large t, so that the energy is conserved. Thus, 

^^i^j + Or 2 ), (18) 
where 5 = r ± t for the initial and final asymptotics, respectively. In the asymptotic 



regions, the spatial Fourier transform of eq.flPSD reads 

^(k) = ^0! / d5(e^ kwkS - e± ikt - ikS )F i>f (8) , (19) 

C i,f 

where k = \k\. The integration contour is ImS = for the final state and lm.5 = T 
for the initial state. 

Consider first the final state. Expanding in e, we find 

where 5 = r — t. The spatial Fourier transform of this field at t — > oo equals 



<P(k) = i_L=.{ e -^*6 k + &I k e^} 

where 



2V37rei? 

= o k 



These amplitudes define the final coherent state, exp(J <ifc6 k 6 k )|0), which has the 
number of particles 



n f = J dkblb k - ' " 



and the energy 



A 



12tt 2 



E = Jdkuo ]l btb ]l = ^-e 2 . (20) 

At t — > — oo, we can find the function Fi(5) in the limit e —>■ 0. In order to find 
the initial state, we need the amplitudes / k (see eq.(f|)). As follows from eq.(|T9D, this 
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coefficient is determined by that part of the function Fi(5) which has singularities 
above the contour ImS = T. The lowest in e contribution of this type is 



16 (5/R - z + ie 2 /8) 2 w y 4 (5/ R - i)(5/ R - i + it 2 /8) 
+ terms, regular at Im 5 > T , 

where 5 = r + t. Calculating the spatial Fourier transform by means of eq.(|19D we 
find 

8^ I / (« + e 2 /8) 2 J 

Eqs.©, (D and (0) give 

where terms subleading in e are omitted. Making use of these equations we obtain 

m= [ dkata k = In 2 / 3 - . (21) 

J A 6 

Note that at small e the number of initial particles is much smaller than the number 
of the final ones. 

The probability of the transition described by the solution fll6|) is given by eq.(|lO|). 
To calculate the imaginary part of the action we note that the only contribution 
comes from the singularity lying between the integration contour and the real time 
axis. Thus, the time integral is reduced to the residue at the pole. The remaining 
space integral can be done exactly. After some algebra one finds 

ImS = S inst = 167T 2 /A 



and, therefore, 

327T 2 

A 

At small e the probability of the transition is exponentially suppressed by the same 
factor as at zero energy. 



- = exp(-^ + 0(e 2 )). (22) 



3 SU(2) model 



In the case of the pure SU(2) theory our strategy of obtaining the solutions is the 
same as in sect. 2. First, we write the spherically symmetric ansatz 

g y 
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where if^ is the t'Hooft symbol and £ = 1/2 ln(?/ 2 / R 2 ). The Yang-Mills equations, 
in terms of the function /(£), read: 



d 2 f 



+ A(f 2 - f)(2f - 1) = 



The solution with the turning point at £ = has the form 




where k = yl — e 2 is a free parameter. In the limit e — > the usual instanton solution 
is reproduced. After the conformal transformation ( |T5D and analytic continuation to 
the Minkowskian domain we obtain[| 

SRtX a f(0 



Ad 
^0 



g[(t - iR) 2 - r 2 }[(t + iR) 2 - r 2 ] 



A a = 4R(S ai (R 2 + t 2 - r 2 ) + 2Re aij Xj + 2x a Xj) 

g[(t - iR) 2 - r 2 ][(t + iR) 2 - r 2 } /lU ' 1 ' 

where r = |jc|. This field configuration clearly has the turning point at t — 0. The 
structure of singularities of this solution is similar to the <p A case discussed in sect. 2. 
In particular, the position of the singularity with smallest Imt is given by eq. fll7D with 
n = 0. 

At first sight, the field seems to violate the asymptotic condition: it does not 
decay as \t\ — > oo. This is, however, merely a gauge artifact (one can check that F^ u 
decays like 1/t). Indeed, the part of the solution that does not go to zero at 
\t\ — > oo is concentrated along the light cone and has the asymptotic form 

x° 2R (1 8-iR\ 



■ 



x S 2 + R 2J \2 5 + iR / 

a _x a x i 2R (I 5 - %R \ 
{ ~ x 2 5 2 + R 2t \2 n 6Tm) ' 

where 8 — r ± t depending on whether the initial or final state is concerned. This 
non-decaying asymptotics can be killed by the gauge transformation with the gauge 
function 

uj = cos F(5) + i<7 a n a sin F(8) , 
1 In fact, this solution was previously obtained in refs.O, E3. 
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where a a are Pauli matrices, n a = x a /r and F(x) satisfies the equation 

dF 2R fl 5-iR\ 

~dJ ~ 5 2 + r 2 * [2 n JTm) ' 

In the new gauge the asymptotics of the transverse spatial components of the field 
23[), which we concentrate on in what follows, reads 



+-{— sin2F + (sin2F - | cos 2F) ^ ^ }{8 ai - n a n t ) 
gr 1 2 R o z + R z } 

Here = - In 5 ~ iR 

In the final asymptotic region we find, expanding in e, 



, , Re 2 ( R 2 (6 2 -R 2 ) 2R 3 S . 

±iW ~ 4 5 x I (5 2 + i? 2 ) 2 + (6 2 + R 2 ) 2 { ai nani) 



where 5 = r — t. Performing the spatial Fourier transform and comparing with the 
plane wave decomposition 



A a ±i (k,t) = 1 ' \eT(k)bl m e- tuJkt + ef{-k)b*^e lu>kt \ 



where e™(&) are the transverse polarization vectors, we obtain 



eT(k)br = eT(-k)b*_T = ^^R 2 e' kR {^e atJ ^ + 5 ai - ^} . (24) 

These expressions can be used to calculate the number particles and the energy of 
the final state. We find 



3e 4 vr 2 



8g 



2 



Consider now the initial state. As in the case of the 4 model, we need to extract 
only that part of the solution which has singularities above the contour Im<5 = T. 
The lowest in e transverse contribution of this type is 

A ^ X) = g^l024 [5/R- i + i e 2 /8W H * ' /R " " ^ " ^ 
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+ terms regular at ImS > T . 
Performing the Fourier transform we obtain 



^-^-^e-^, (26) 



'QAgVk [ mJ k k 2 



where the last estimate is obtained at k <C 1/e 2 . From eqs.([F[), ® and ( p5|) we find 

T - AT = R\l - f^) 1/3 \ . 



Making use of this expression we calculate the total number of initial particles, 

As in the case of the 4 model, at small e the number of particles in the initial 
state is much smaller than the number of particles in the final state. Note also that 
the characteristic momentum of the initial particles is k in ~ e~ 4 / 3 <C e~ 2 , so the 
assumption of eq.(^) is justified. 

Our next step is to calculate the imaginary part of the action. As in the case of 
the 4 model, we first evaluate the integral over t along the contour of fig. 2. The 
only contribution to the imaginary part comes from the pole which lies between the 
contour AB and the real time axis. The remaining integral over d 3 x can be done 
exactly and gives 

ImS* 
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T 

i.e. the usual instanton action. Estimation of other contributions into eq. (PD finally 
gives 

a = exp{-^ + 0(e 2 )} (27) 

Again, at small e the probability of the transition is exponentially suppressed by the 
same factor as at zero energy. 

Unlike the case of 4 model, we can explicitly check that our solution interpo- 
lates between two different topological sectors. This can be done by calculating the 
topological charge along the contour ABCD of fig.l. Since the solution is symmetric 
with respect to t — > —t, the topological charge evaluated along the real time axis is 
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zero. Therefore, the only contribution to the topological charge comes from the pole 
lying between the contour AB and the real time axis. However, as one can explicitly 
check, for our solution the residues of (F^F^) and {—F^ V F^ U ) at this pole coincide, 
i.e. the solution is anti-self-dual at the singularity. Thus, the topological charge is 

Q = -f^lmS = -1 , 
and our solution indeed interpolates between different topological sectors. 

4 Conclusion 

To summarize the results of sects. 2 and 3, we have obtained the set of solutions to the 
classical field equations that satisfy the requirements i) and ii) and, therefore, saturate 
the probability of the instanton-like transition from some coherent state \a) in the 
saddle point approximation. In both (f) 4 model and the SU(2) gauge theory the set of 
solutions depends on two parameters, the scale parameter R and the dimensionless 
parameter e. In the limit of small e the initial coherent state and the most probable 
final state for the corresponding process were found. In both models the number of 
initial particles is much less than the number of final particles, while the latter is much 
smaller than the characteristic number of particles in the sphaleron configuration (1/A 
and l/g 2 ). The characteristic momenta of particles in the final state are determined 
by the scale parameter, kf ~ 1/R, while the characteristic momentum in the initial 
state is much larger, ki ~ e _4//3 /i?, in both models. Thus, in the limit e — > any of 
our solutions describes an instanton-like transition from a particular coherent state 
containing small number of particles at the energy given by eqs.(^) and fl2"5|). The 
probability of this transition goes to exp(— 2S inst ) at e — ► 0. 

One may ask whether the solutions ( |TED and ( |2"3"D saturate also some transition 
probability maximized over the initial states, i.e. the quantity J2ab I (&|»S'-Pe-Pa|o) | 2 , 
where Pa is the projector onto the eigenspace corresponding to the eigenvalue A of the 
operator A = / dk\(k)a\ i a\ s ,. The criterion is |TD[ that the ratio /k//k of the positive- 
and negative-frequency amplitudes in the initial state is const • exp(— A^c) where c is 
some real number. Our solutions fll6|) and ( Effi) do not satisfy this condition at any 
real Ak- The reason is that for both 4 model and SU(2) gauge theory the product 
/k/k changes the sign at some k, in exact analogy to the case of 0(3) sigma-model 
[§]. Therefore, our solutions cannot be interpreted as saturating some maximum 
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probability and does not automatically set an upper bound to the two-particle cross 
section. 

To draw some physical conclusions from the existence and properties of the solu- 
tions and (^), we have to assume that i) the probability of the transition does 
not depend, with exponential accuracy, on the details of the initial state in the limit 
of small number of initial particles and ii) there are no other solutions describing the 
same process with greater action. If these assumptions were justified, the two-particle 
cross section at high energies would coincide (with exponential accuracy) with the 
probability of the transitions described by the solutions flT6| ) and ( p3|) in the e — > 
limit and thus would be suppressed at high energies by the same suppression factor 
exp(— 2Si ns t) as at zero energies. 

As discussed in the Introduction, the first assumption is supported by the per- 
turbative calculations at low energies. Whether the second assumption is justified 
or not is presently unclear. This question is closely related to the (equally unclear) 
following one: at what energies the massless theories discussed in sects. 2 and 3 can 
be considered as an approximations to the massive ones? Naively, one would expect 
that the massless approximation is valid at least at mR <C 1 since R is the largest 
scale. In view of eqs.(|20|) and fl25|) this leads to the conditions 

E ^> e 2 E sp h in the 4 model , 

(28) 

E > e A E sph in the SU(2) theory , 

instead of the expected E ^> E sp h- 

If the conditions fl28|) are correct, they cause serious problems with the interpre- 
tation of our solutions. Since these conditions overlap with the region of validity of 
the perturbation theory around the zero energy instanton, it is a crucial point to see 
whether our results, eqs. (p2|) and (0), can be reproduced by perturbative methods 
when we take as an input the initial states of sects. 2 and 3. The perturbative saddle 
point at small number of initial particles was studied in ref.[0. It was found that, in 
the leading order approximation, the probability of the transition in that case grows 
exponentially with energy and thus differs from the behaviour we have found in sects. 2 
and 3. Simple estimates show that for our initial states there exists another saddle 
point (in addition to that described in ref.0) with the action and number of particles 
which parametrically coincide with those found in sects. 2 and 3. If this estimate was 
reliable, we would conclude that our saddle point is irrelevant at least at low energies, 
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since it has smaller action than the saddle point of ref . |7j] . However, the perturbative 
expansion about the new saddle point is ill defined, so that the saddle point itself 
may be an artifact of the leading order approximation. Thus, we do not arrive at 
the self-consistent picture and cannot make definite conclusions about the relevance 
of our solutions to the instanton-like transitions at the moment. 

The authors are indebted to V.A.Rubakov for numerous and helpful discussions. 
The work of D.T.S. is supported in part by the Weingart Foundation through a 
cooperative agreement with the Department of Physics at UCLA. 
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FIGURE CAPTIONS 



Fig.l The contour in the complex time plane where the solution corresponding to the 
instanton-like transition is defined. A and D mark initial and final asymptotic 
regions, respectively. The line BC represents the Euclidean part of the evolution. 

Fig. 2 The structure of singularities of the solution in the complex time plane. Solid 
curved lines represent the positions of the poles as \x\ chagnes from oo to 0. 
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